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Abstract 

Options are financial instruments that depend on the underlying 
stock. We explain their non-Gaussian fluctuations using the nonexten- 
sive thermodynamics parameter q. A generalized form of the Black- 
Scholes (B-S) partial differential equation, and some closed-form so- 
lutions are obtained. The standard B-S equation (q = 1) which is 
used by economists to calculate option prices requires multiple values 
of the stock volatility (known as the volatility smile). Using q = 1.5 
which well models the empirical distribution of returns, we get a good 
description of option prices using a single volatility. 

Although empirical stock price returns clearly do not follow the log- 
normal distribution, many of the most famous results of mathematical fi- 
nance are based on that distribution. For example, Black and Scholes (B-S) 
[1] were able to derive the prices of options and other derivatives of the un- 
derlying stock based on such a model. An option is the right to buy or sell 
the underlying stock at some set price (called the strike) at some time in the 
future. While of great importance and widely used, such theoretical option 
prices do not quite match the observed ones. In particular, the B-S model 
underestimates the prices of options in situations when the stock price at the 
time of exercise is different from the strike. In order to match the observed 
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market values, the B-S model would need to use a different value of the 
volatility for each value of the strike. Such "implied volatilities" of options of 
various strike prices form a convex function known as the "volatility smile" . 

Indeed, attempts have been made to modify the B-S model in ways that 
can correct for the smile effect (cf [2] or more recently [3, 4]). However, those 
approaches are often very complicated or rather ad-hoc, and do not result 
in managable closed form solutions, which is the forte of the B-S approach. 
In this paper we do however succeed in developing a theory of non-Gaussian 
option pricing which allows for closed form solutions for European options, 
which are such that can be exercised exclusively on a fixed day of expiration 
and not before (as is the case for American options). 

Our approach uses stochastic processes with statistical feedback [5] as a 
model for stock prices. Such processes were recently developed within the 
Tsallis generalized thermostatistics [6]. The driving noise can be interpreted 
as a generalized Wiener process governed by a Tsallis distribution of entropic 
index q. In the limit q — > 1 the standard model is recovered. For q « 1.5, 
this model closely fits the empirically observed distribution for many financial 
time series, such as stock prices [7] (Figure 1), SP500 index, [7]) [8], FX rates, 
etc. This is consistent with a cumulutive distribution having power tails 
of index 3 [9]. We derive closed form option pricing formulas, reproducing 
prices which, relative to the standard B-S model, exhibit volatility smiles 
very close to those observed empirically (Figure 4). Note that q — 1.5 well 
models hydrodynamic turbulence on small scales [10], reinforcing notions of 
a possible analogy between these two systems. 

The standard model for stock prices is S = S e Y ^ where Y(t) = ln(S(t + 
to)/ 4n S(t ) follows 

dY = fidt + adu (1) 

The drift fi is the mean rate of return and a 2 is the variance of the stock 
logarithmic return. The noise w is a Brownian motion defined with re- 
spect to a probability measure F. It is a Wiener process and satisfies 
E F [du(t)duj(t')} = dtdt'S(t — t') where the notation E F \\ means the expec- 
tation value with respect to the measure F. This model yields a Gaussian 
distribution for Y resulting in a log-normal distribution for S. Within this 
framework, Black and Scholes were able to establish a pricing model to obtain 
the fair value of options on the underlying stock S. 

In this paper we assume that the log returns Y(t) = \nS(t + t )/\nS(t ) 
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follow 

dY = fxdt + adtt (2) 

with respect to the timescale t. Here Q evolves according to the statistical 
feedback process[5] 

dQ = Pin^du (3) 
The probability distribution P satisfies the nonlinear Fokker-Planck equation 

|p(n,t|n',o = ^P 2 - ff (n,*|n / ,0 (4) 

Explicit solutions for P are given by Tsallis distributions [11] 

p q (n,t | n(o),o) = ^{1 - - q)Mt) - m?}^ (5) 

Choosing (5{t) = c^{(2-q)(3-q)t)- 2 ^ 3 -^ and Z(t) = ((2-g)(3-g)ct)s=i 
ensures that the initial condition P q = S(Q(t) — fi(0)) is satisfied. The q- 

dependent constant c is given by c = f3Z 2 with Z = j ao 0O {\ — {\—q)l3VL 2 )~dVt 
for any (3. With £1(0) = 0, we obtain a generalized Wiener process, dis- 
tributed according to a zero-mean Tsallis distribution In the limit q — > 1 
the standard theory Eq(l) is recovered, and P q becomes a Gaussian. We 
are concerned with the range l<g<5/3in which positive tails and finite 
variances are found [12] . The distribution for In S becomes 

P q (\nS(t+t ),t+t | \nS(t ),t ) = ^{l-^(t)(l-g)[ln ^+j o) -/it] 2 }^ 

(6) 

with (3 = (3(t)/a 2 . This implies that log-returns \n[S(t + t )/S(t )] over the 
timescale t follow a Tsallis distribution, consistent with empirical evidence 
for several markets, e.g. S&P500 (Figure 1 [7]) [8], with q « 1.5. 

Our model exhibits a feedback from the macroscopic level characterised 
by P, to the microscopic level characterised by Q. We can imagine that this 
is really due to the interactions of many individual traders whose actions all 
contribute to shocks to the stock price which keep it in equilibrium. Their 
collective behaviour yields a nonhomogenous reaction to returns: rare events 
(i.e. extreme returns) will be accompanied by large reactions, and will tend 
to be followed by large returns in either direction. 
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Using Ito calculus [13, 14], the equation for S follows from Eq(2) as 



dS = jlSdt + aSdQ (7) 

where jl = /i+^-P^ 1-9 . The term £ ^-P q l ~ q is the noise induced drift. Remember 
that P q is a function of Q with 

n«) = l " s(f)/ll '/ (0) -' ri (8) 

(with to = for simplicity.) As in the standard case (cf [2]), the noise term 
driving S is the same as that driving the price f(S) of a derivative of the 
underlying stock. It should be possible to invest one's wealth in a portfolio of 
shares and derivatives such that the noise terms cancel each other, yielding 
a risk- free portfolio, the return on which must be the risk-free rate r. This 
results in a generalized B-S PDE 

^L + rS^ + -^ a 2 S 2 P l - q = rf (9) 
dt + dS + 2d£P q 1 [) 

where P q (Q(t)) evolves according to Eq(4). For q — > 1 the standard case 
is recovered. This PDE depends explicitly only on the risk-free rate and 
the variance, not on //, but it does depend implicitly on /i through P q (£l), 
with given by Eq(8). Therefore, to be consistent with risk-free pricing 
theory, we should first transform our original stochastic equation for S into 
a martingale before we apply the above analysis. This will not affect our 
results other than that jl will be replaced by the risk-free rate r, ultimately 
eliminating the dependency on /i. We now show how this is done. 

The discounted stock price G = e~ rt S follows dG = (jl — r)Gdt + aGdVL 
where dVL follows Eq(3). For there to be no arbitrage opportunities, risk-free 
asset pricing theory requires that this process be a martingale, which it is 
not due to the drift term (jl — r)Gdt. One can however define an alternative 
driving noise z associated with an equivalent probability measure Q so that, 
with respect to the new noise measure, the discounted stock price has zero 
drift and is thereby a martingale. Explicitly, 

dG = (jl-r)Gdt + aGP^du (10) 

Here, P is a non-vanishing bounded function of Q. With respect to the initial 
noise u>, f2 relates to S via Eq(8). That is why for all means and purposes, 
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P in Eq(10) is simply a function of S (or G), and the stochastic process 
can be seen as a standard state-dependent Brownian one. As a consequence, 
both the Girsanov theorem (which specifies the conditions under which we 
can transform from the measure F to Q) and the Radon-Nikodym theorem 
(which relates the measure F to Q) are valid, and we can formulate equivalent 
martingale measures much as in the standard case [15, 16, 17]. We rewrite 
Eq(10) as 

dG = oGP^dz (11) 
where the new driving noise term z is related to u through 

dz = ^~j} dt + dw (12) 
aP q 2 

With respect to z, we thus obtain dG = aGdVt with dVt = P q 2 dz which is 
non other than a zero-mean Tsallis distributed generalized Wiener process, 
completely analogous to the one defined in Eq(3). Transforming back to S 
yields dS = rdt + aSdQ. Compared with Eq(7), the rate of return ft, has been 
replaced with the risk-free rate r. This recovers the same result as in the 
standard asset pricing theory. Consequently, in the risk-free representation, 
Eq(8) becomes 

^S-'^iVw)^) (13) 

This eliminates the dependency on fj, which we alluded to in the discussion 
of Eq(9). As discussed later on, by standardizing the distributions P q (Q(s)) 
we can explicitly solve for Q(t) as a function of S(t) and r. 

Suppose that we have a European option C which depends on S(t), whose 
price / is given by its expectation value in a risk-free (martingale) world as 
f(C) = E® [e~ rT C] . We assume the payoff on this option depends on the 
stock price at time T so that C = h(S(T)). After stochastic integration of 
Eq(ll) to obtain S(T) we get 

/ = e ~ rT E Q 
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h(s{Q)e^U\p q l7 dz s + j\r-^P^)ds\) (14) 



1-q 

The key point is that the random variable J P q 2 dz s = J dQ(s) = fi(T) 
follows the Tsallis distribution Eq(5). This gives 



—rT 

e 



S(0) exp(afi(T) + rT - —aT*k + (1 - q ) a T^^^a 2 n 2 {T)) 



(l-P(T)(l-q)n(Tf) — dn T (15) 

with a = |(3 — q)((2 — q)(3 — q))c) 3 -i . We have utilized the fact that each 
of the distributions P(Q(s)) occuring in the latter term of Eq(14) can be 
mapped onto the distribution of fi(T) at time T via the appropriate variable 



transformations Q(s) = yJ(3(T)/(3(s)Q(T). A major difference to the stan- 
dard case is the f2 2 (T)-term which is a result of the noise induced drift. With 
q — 1, the standard option price is recovered [15]. 

Eq(15) is valid for an arbitrary payoff h. We shall evaluate it explicitly 
for a European call option, which gives the holder the right to buy the stock 
S at the strike price K, on the day of expiration T. The payoff is C = 
m&x[S(T) — K, 0]. Only if S(T) > K will the option have value at expiration 
T (it will be in-the-money). The price c of such an option becomes 

c = E Q [e~ rT C} = E Q [e- rT S(T)] D - E Q [e~ rT K\ D = J x - J 2 (16) 

where the subscript D stands for the set {S(T) > K}. This condition is 
met if -^-aT^i + (1 - q)aT^^a 2 n 2 + atl + rT > \nK/S(0), which is 
satisfied for Q between the two roots si and s 2 of the corresponding quadratic 
equation. This is a very different situation from the standard case, where the 
inequality is linear and the condition S(T) > K is satisfied for all values of 
the random variable greater than a threshold. In our case, due to the noise 
induced drift, values of S(T) in the risk- neutral world are not monotonically 
increasing as a function of the noise. As q — > 1, the larger root s 2 goes toward 
oo, recovering the standard case. But as q increases, the tails of the noise 
distribution get larger, as does the noise induced drift which tends to pull 
the system back. As a result we obtain 

Ji = ^(0)^y£ 2 exp(^-^ a T^-(l-g)aTA^a 2 fi 2 ) 

(1 - (1 -q)(3{T)Q 2 )^dQ (17) 
J 2 = e- rT K^J s S \l-(l-q)(3(T)n 2 )^dn (18) 
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The equation Eq(16) with Eq(17) and Eq(18) constititutes a closed form 
expression for the price of a European call. We calculated option prices for 
different values of the index q, and studied their properties as a function of 
the relevant variables such as the current stock price S'(O), the strike price 
K, time to expiration T, the risk free rate r and o. The results obtained 
by our closed form pricing formula were confirmed both by implicitly solving 
the generalized B-S PDE Eq(9) and via Monte Carlo simulations of Eq(ll). 
Note that American option prices can be solved numerically via Eq(9). 

We compare results of the standard model (q = 1) with those obtained 
for q = 1.5, which fits well to real stock returns. Figure 2 shows the differ- 
ence in call price. In Figure 3, the B-S implied volatilities (which make the 
q = 1 model match the q = 1.5 results) are plotted as a function of K. The 
assymetric smile shape, which is more pronounced for shorter times, repro- 
duces well-known systematic features of the "volatility smile" that appears 
when using the standard q = 1 model to price real options. In Figure 4, 
the volatility smiles for actual traded options on BP and S&P 500 futures is 
shown together with those resulting from our model using q = 1.5. These 
results are encouraging, and we are currently studying a larger sample of op- 
tions data. Empirical work is required to see if arbitrage opportunities can 
be uncovered that do not appear when the standard model is used. Another 
potential application will be with respect to option replication and hedging. 
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10 NASDAQ stocks, 1 min, 2001 

10° 1 1 1 1 1 




Normalized Return 



Figure 1: Distributions of log returns over 1 minute intervals for 10 high- volume stocks, 
normalized by the sample standard deviation. Also shown is the Tsallis distribution of 
index q = 1.43 (solid line) which provides a good fit to the data (Figure kindly provided 
by R. Osorio). 



9 



Call Price Difference 




' 20 30 40 50 60 70 80 

Strike K 

Figure 2: Calibrated so that the options are priced equally for 5(0) = K, the difference 
between the q = 1.5 model and the standard B-S model is shown, for S(0) — $50 and 
r = 0.06. Solid line: T = 0.6 with a = .3 for q = 1 and a = .297 for q = 1.5. Dashed line: 
T = 0.05 with a = .3 for q = 1 and a = .41 for q = 1.5. Times are expressed in years, r 
and a are in annual units. 
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Implied Volatility 




option prices, one can back out the volatilities implied by a standard q = 1 B-S model. 
T = 0.1 (circles), T = 0.4 (triangles). 
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a)lmplied Vol(BP) 



b)lmplied Vol(SX) 




Strike K Strike K 



Figure 4: A comparison of option prices from a q = 1.5 model and traded prices is 
given by a comparison of volatility smiles. a)Implicd vols, for options on British Pound 
futures (last trading date Dec 7, 2001) [18] vs. strike (for 5(0) = 143, r = .065, T = .0055 
(2 days)) (symbols); implied vols needed for a q — 1 B-S model to match prices from a 
q = 1.5 model using a = .1445 across all strikes (line).b) Implied vols for S&P500 futures 
(SX June, last trading date June 15, 2001, 5(0) = 1275, r = .065, T = .027 (10 days)) 
(symbols); implied vols, from q = 1.5 model with a — .3295 (line). 
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